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Abstract 

The final state of turbulent magnetic relaxation in a reversed field pinch is well explained by 
Taylor's hypothesis. However, recent resistive-magnetohydrodynamic simulations of the relaxation 
<^ \ of braided solar coronal loops have led to relaxed fields far from the Taylor state, despite the 

conservation of helicity. We point out the existence of an additional topological invariant in any 
flux tube with non-zero field: the topological degree of the field line mapping. We conjecture that 
this constrains the relaxation, explaining why only one of three example simulations reaches the 
Taylor state. 
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PACS numbers: 52.30.Cv, 52.65.Kj, 52.35.Vd, 96.60.Hv 



The landmark paper of J. B. Taylor |l| showed how the final, relaxed, magnetic field in 
a turbulent plasma experiment could be predicted theoretically by assuming conservation 
of the total magnetic helicity H = J v A • B d 3 x, where B = V x A. Taylor hypothesized 
that, in a turbulent plasma with small but no n- vanishing resistivity, H is the only relevant 
constraint on the relaxation. The resulting minimum energy state is then a linear force-free 
field [2|, V x B = <xB, where a is a constant depending on the value of H. 

The success of this theory in predicting the final state of relaxation in a reversed field pinch 
has led to speculation that Taylor's hypothesis might apply more generally, for example, to 
magnetic structures in the Sun's atmosphere J3|-l6|. Here, the nature of the final state is of 
great interest because it limits the magnetic energy available for conversion to heat during 
the relaxation. This is key to understanding how magnetic fields produce the extreme coronal 
temperatures in the Sun and other stars. To apply Taylor's hypothesis in such magnetically 
open domains, H is replaced by the relative helicity Hr with respect to a reference field 
|4, [Z|, typically a potential field. A local domain, such as an isolated loop, is taken (the 
corona cannot be globally a linear force-free field) and the hypothesis has met with some 
notable success |6j. However, counter-examples are known |8|, |9j and it now seems likely 
that constraints beyond the conservation of global magnetic helicity may be important for 
general relaxation events. For example, while Hr measures only the second-order linkage 



of field lines, higher-order invariants may play a role |l0l-ll3j (although this remains under 
debate \l£ 



This Letter presents a conjecture on one possible topological constraint and was motivated 
by recent resistive-magnetohydrodynamic simulations of the relaxation of a magnetic flux 
tube, intended to model a solar coronal loop. In all of our experiments the initial flux tube 
contains a braided magnetic field of non-trivial topology but with no net current or helicity 
[9|, |l5[ . Despite undergoing a turbulent relaxation with myriad small-scale current structures 
and magnetic reconnections, on a timescale short compared to the resistive timescale, the 
magnetic field relaxes to a final state which directly contradicts that predicted by Taylor's 
hypothesis. Although force- free, the final state is not linear force- free, because the coefficient 
a varies strongly across field lines. For the first experiment, where we started from a field 
modeled on the pigtail braid, the current is concentrated into two flux tubes of opposite 
twist (sign of the parallel current). The initial and final states of this braiding experiment 
are shown in Fig. [TJ If Hr were the only relevant constraint on relaxation then the final state 
would be a near- uniform, vertical, potential field (B = V"0), because the initial configuration 
is chosen such that Hr = 0. Note that Hr is still well conserved during the turbulent 
relaxation, so the departure of the final state from the Taylor prediction is not due to 
changes in Hr. Rather, the result seems to indicate the presence of additional constraint (s). 

In this Letter, we propose that such a constraint is given by a property of the field line 
mapping called its topological degree. This is a topological property that can be defined for 
any continuous mapping from one compact manifold to another. In the braiding experiment, 
the mapping of magnetic field lines from the lower boundary So (where they enter the 
domain) to the upper boundary Si (where they exit) forms such a mapping. This is because 
(a) all field lines connect, in the same direction, these two boundaries (the vertical flux 
through any horizontal cross-section is constant through the domain), and (b) the magnetic 
field is everywhere non-zero (B 7^ 0) within the domain. Although the fields considered in 
this paper are defined in rectangular, Cartesian domains such that the flux tube is aligned 
along the z-axis and the boundaries are rectangles in the (x, y) plane, the ideas are more 
generally applicable because every flux tube of any shape with uni-directional field B/0 
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FIG. 1. (Color online) Magnetic field lines in the original braiding experiment 15j for (a) the 
initial state at t = 0, and (b) the relaxed state at t = 290 (in units of the AlfVen time). Field lines 



are traced from the same starting points in each case. In (b), color contours show a 
the Sq plane. 
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is topologically equivalent to a straight cylinder [161 ] . 

To define the topological degree of the field line mapping, let F(x ) G Si be the end-point 
of the field line starting at the point x 6 <Sn. The field line mapping will have one or more 



periodic orbits Xq G S where F(xq) = x£ 



a: 



As a fixed point of the mapping F — I (where 



I denotes the identity map), each Xg is characterized by a fixed point index, denned as the 
local Brouwer degree of the mapping (see e.g. [la]). This index takes integer values, ±1 for 
generic, structurally stable, isolated periodic orbits. The case +1 corresponds to an elliptic 
null point of the local linearisation of F — I, or an elliptic periodic orbit, while —1 indicates 
a hyperbolic periodic orbit. The sum 



T = y^ index (xq) 



(1) 



over all isolated periodic orbits is called the Lefschetz number or topological degree of the 



mapping F. By the Lefschetz-Hopf theorem |19f] it is a conserved quantity, providing that 



no periodic orbits cross the side boundary of the domain. Periodic orbits can therefore be 
created or annihilated only in pairs of opposite index. 

The topological degree T may be computed by evalua ting the Kronecker integral around 
the boundary of S [20J, or by other numerical methods 2l| . Here, we adopt the graphical 



t= 0. 



t= 35. 




FIG. 2. (Color online) Sequence of color maps for the original braiding experiment 15(, where 
total degree T = 2. 



color map technique of [20] . For example, Fig. [2] shows the colour maps at various times 
in the braiding simulation. Every point xo = (xq, yo) on the lower boundary of the domain 
(So) is assigned one of four colours, according to its field line mapping F(x ) = (F x , F y ) to 
the upper boundary Si. We use red if F x > xq and F y > y ; yellow if F x < xq and F y > y ; 
green if F x < x and F y < y ; and blue if F x > x and F y < y . On the resulting color map, 
periodic orbits correspond either to red-green or yellow-blue boundaries; isolated, generic 
periodic orbits are points where all four colours meet. The index of an isolated periodic 
orbit may be read from the sequence of colors passed through in an anti-clockwise direction 
around a small circle around the point. An elliptic orbit (index 1) has red-yellow-green-blue 
(r-y-g-b), while a hyperbolic orbit (index —1) has r-b-g-y. To determine T, either sum the 
individual indices or simply record the sequence of colors around the boundary of So- For 
example, in Fig. E^a), the initial state for the braiding experiment, we find 12 periodic orbits 
with index 1, and 10 with index -1, giving a net topological degree of 2. This corresponds 
with the (anti-clockwise) sequence of colors on the boundary of r-y-g-b-r-y-g-b . 

The later color maps in Fig. [2] for the braiding experiment show that, during the turbu- 
lent relaxation, changes in the magnetic topology by reconnection lead to the annihilation 
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FIG. 3. (Color online) Color map for the initial magnetic field of the Browning et al. [y] simulation 
(case 1). 

of periodic orbits in pairs of index +1 and —1. This is consistent with (JT|) as the reconnec- 
tion processes occur strictly in the interior of the domain such that T remains unchanged. 
Eventually only two elliptic periodic orbits remain in the final relaxed state, lying within 
the two flux tubes of opposite twist shown in Fig. []Jb). This leads us to suggest an expla- 
nation for the discrepancy between this relaxed state and that predicted by Taylor theory: 
for a turbulent relaxation in a magnetic flux tube that leaves the side boundary fixed, the 
topological degree T of the initial field line mapping is an additional topological constraint 
that can prevent the system from reaching the Taylor state. 

We now consider another two resistive-MHD simulations of turbulent relaxation in mag- 
netic flux tubes which further support our degree conjecture. The first is a simulation by 
Browning et al. [6J where the Taylor hypothesis succeeds in predicting the relaxed state. 
Using their expressions for B, we have computed the color map (shown in Fig. [3]) of the 
initial force-free magnetic field in their case 1. The field consists of an axisymmetric twisted 
magnetic flux tube of unit radius, with a piecewise-constant profile of a with radius. Fig. [3] 
shows that this initial configuration has a single isolated periodic orbit at the origin. In 
addition, there are four rational surfaces on which field lines have a winding number that 
exactly divides 2-k, forming four continuous circles of periodic orbits in the So plane. How- 
ever, it is easy to show that such rational surfaces do not contribute to the total degree, so 
T = 1, in accordance with the pattern of colors around the tube boundary, r = 1. 

The initial flux tube in the Browning et al. simulation is kink unstable, and a turbulent 
relaxation is triggered by an initial velocity perturbation. The system relaxes to a force-free 
equilibrium on a comparable timescale to our original braiding experiment. However, the 
difference here is that this final equilibrium is a constant-a "Lundquist" magnetic field, in 
accordance with expectations from Taylor theory |22|. The simulation still supports our 
degree conjecture because the final state retains a total degree T — 1, with a single twisted 
flux tube. The difference between this and our braiding experiment is that here the Taylor 
state is compatible with the degree of the initial field. 

Our final example is a new simulation of an initially braided magnetic field with degree 



T = 3, solving the same resistive-MHD equations as in our earlier braiding simulation [15 



The orig inal initial condition with T = 2 was modeled on the pigtail braid (see Fig. 2 of 



Ref. (23J), and was formed from a uniform vertical field super-imposed with six localized 



"twist" regions (Eq. 2 of [23j) arranged in two columns. The T = 2 property arose from 



the arrangement of the positive twist regions in one column and the negative regions in 
the other. We obtain a topologically stable T = 3 configuration by using four columns of 
twist regions arranged on a circle, where each contribute an elliptic periodic orbit. The total 
degree is 3 because the remainder of the field contributes a net degree of —1. The initial field 
we used is shown in Fig. Ufa). In the notation of Ref. [23J, the twist regions have param- 
eters ci = (r , 0,-18, 1,^2, 2), c 2 = (-r , 0, -18, 1,^2,2), c 3 = (0,^,-6,-1,^,2), 
c 4 = (0,-7-0,-6,-1,^,2), c 5 = (r , 0,6, 1,^2, 2), c 6 = (-r ,0, 6, 1, ^2,2) , c 7 = 
(0,r , 18, -1,^,2), and c 8 = (0,-^,18,-1,^,2), where r = 1.27. Note that the 
total helicity of this field vanishes, as in the T = 2 experiment. The color map for this 
initial field (Fig. Hb) confirms that T — 3. 

The T = 3 initial condition leads again to a turbulent relaxation, reaching a relaxed state 
on a similar timescale to the T = 2 experiment. The final color map is shown in Fig. H](d). 
As in the T = 2 case, the Taylor theory fails, but the degree conjecture is upheld, with 
T = 3 maintained in the final state. In fact, there are four elliptic periodic orbits and one 
hyperbolic periodic orbit, so the four initial columns have led to four twisted flux tubes 
in the final state. The degree constraint itself would not preclude the annihilation of the 
hyperbolic orbit with one of the elliptic orbits. There must be some further topological 
constraint on the evolution preventing this. This remains to be fully explored. 

In conclusion, our examples show that preservation of the topological degree T leads to 
non-trivial constraints on braided magnetic flux tubes. We therefore claim that the topolog- 
ical degree imposes a constraint on the relaxation beyond that of helicity conservation. In 
particular, the Taylor state is not reached if its degree differs from that of the initial state, 
at least in cases like our examples where the boundary is unaffected by the dynamics. Note 
that, although these simulations use line-tied boundary conditions (in z), the idea extends 
to the case of a ^-periodic boundary and hence to a toroidal domain. Finally, we point out 
that a series of further invariants may be constructed by taking the topological degree of 
multiple iterations of the field line mapping, that is, by considering the degree of periodic 
orbits with periods greater than one. 
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